Expansion of an interacting Fermi gas 
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We study the expansion of a dilute ultracold sample of 
fermions initially trapped in a anisotropic harmonic trap. The 
expansion of the cloud provides valuable information about 
the state of the system and the role of interactions. In par- 
ticular the time evolution of the deformation of the expand- 
ing cloud behaves quite differently depending on whether the 
system is in the normal or in the superfluid phase. For the 
superfluid phase, we predict an inversion of the deformation 
of the sample, similarly to what happens with Bose-Einstein 
condensates. Viceversa, in the normal phase, the inversion of 
the aspect ratio is never achieved, if the mean field interaction 
is attractive and collisions are negligible. 

Since the first experiments on trapped Bose-Einstein 
condensates the imaging of the expanding cloud, fol- 
lowing the sudden switching off of the confining poten- 
tial, has provided crucial information on the novel fea- 
tures exhibited by atomic gases in conditions of quan- 
tum degeneracy. These include, in particular, the bi- 
modal structure of an expanding Bose gas at finite tem- 
perature and the anisotropy of the asymptotic profile of 
the condensate [1,2]. In the Thomas-Fermi limit, where 
the Gross-Pitaevskii equations coincide with the hydro- 
dynamic theory of superfiuids, the anisotropy of the ex- 
panded gas reflects the anisotropy of the pressure force 
which is stronger in the direction of tighter confinement 
[3] . The predictions of the hydrodynamic equations and 
of the consequent scaling behaviour exhibited by Bose- 
Einstein condensates during the expansion have been in- 
vestigated by several authors ( [4], [5], [6]), providing ex- 
cellent agreement with experiments [7,8] and pointing out 
the difference with respect to the expansion of a non con- 
densed gas. In the latter case, in the coUisionless regime, 
the density profile approaches an isotropic shape, inde- 
pendent of the initial deformation of the gas. 

In this letter we study the problem of the expansion 
of an ultracold sample of fermions initially trapped in 
an anisotropic harmonic trap. We will show that also in 
the case of fermions the expansion of the gas provides 
valuable information about the state of the system and 
the role of interactions. We will consider a gas of atoms 
interacting with attractive forces. This is a natural re- 
quirement for the realization of Cooper-pairs and hence 
for the achievement of the superfluid phase [9] . Such in- 
teractions are naturally present in some fermionc species 
like ^Li and can otherwise be obtained by changing the 



scattering length profiting of the presence of a Feshbach 
resonance. 

The description of the expansion of a cold fermionic 
gas in the normal and superfluid phase requires differ- 
ent theoretical approaches. For the normal phase we use 
the formalism of the Landau- Vlasov equations, while in 
the superfluid phase we study the expansion using the 
hydrodynamic theory of superfiuids. 

We consider the case of two different fermionic states, 
hereafter called 1 and 2, initially confined in a harmonic 
trap. We assume that the two species are present in 
the same amount and feel the same trapping potential, 
so that the densities of the two species are equal at 
equilibrium as well as during the expansion: ni(r, i) = 
ri2(r, i) = rt(r,t)/2. The trapping potential will be cho- 
sen of harmonic type 
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describing a cilyndrically simmetric trap with deforma- 
tion A = lJz/^Jl- The interaction between the two 
fermionic species is fixed by the coupling constant g — 
ATTh^a/m, where a is the s-wave scattering length. 
In this paper we will use the equation of state 
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to describe the uniform phase of the gas where the first 
term is the kinetic energy evaluated at zero temperaure 
while the second one is the interaction energy evaluated 
in the mean field approximation. Equation (2) neglects 
the effects of correlations which become important for 
large values of the scattering length and affect in a dif- 
ferent way the equation of state of the normal and su- 
perfluid phase. The formalism developed in this paper 
can be easily generalized to include a more accurate de- 
scription of the equation of state. It is however worth 
pointing out that, even using the same equation of state 
for the normal and superfluid phases, the expansion of 
the gas behaves quite differently in the two cases being 
described by different kinetic equations. In the presence 
of the external potential (1), the equilibrium condition 
in the local density approximation is determined by the 
equation 
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where /i is the chemical potential of the sample fixed by 
the normalization condition. 

The relevant parameter characterizing the interaction 
in the fermionic system is the ratio 



X = 



int 



Eho 



between the interaction energy 

and the oscillator energy 

Eho = J Vho{r)n{r)(fr. 
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In the perturbative regime the integrals (5,6) can be eval- 
uated using the non interacting density profile which, in 
Thomas- Femi approximation (3), takes the simple form 
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After integration of (5) and (6) one finds [10] 
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where kp = (37r^n(0))^/^ is the Fermi momentum eval- 
uated at the central value of the density. In order to go 
beyond the perturbative regime, one determines numeri- 
cally the ground state density by minimizing the energy 
of the system 
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+ J Vho{r)n{T)d^r + ^ J n\v)d'r, (9) 



One easily finds that the equilibrium value of the ra- 
tio Eint /Eho depends on the dimcnsionlcss combination 
aN^/^/aho also in the non perturbative regime [9], as 
reported in Fig.l. Using (9) one predicts that the com- 
pressibility of the gas becomes negative in the center of 
the trap if \a\N^/^ /aho > 0.61. For large values of \a\ the 
resulting predictions should be however taken with care 
since the functional (9) ignores correlation effects beyond 
mean field. 

Let us now discuss the expansion of a fermionic sam- 
ple trapped in an elongated harmonic trap (A < 1). We 
describe first the expansion of the normal fluid and after- 
wards the one of the superfluid. In the ideal case, using 
the semiclassical description, one finds that the ratio of 
the square radii evolves according to the classical law 
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The ratio (10) approaches unity for large times, reflect- 
ing the isotropy of the momentum distribution. This 
result ignores the effects of collisions which are however 
expected to play a minor role at low temperature, due to 
Pauli blocking, unless the scattering length is very large. 
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FIG. 1. Ratio x = Eint/ Eho as a function of the universal 
parameter aN^^^ / auo calculated with the mean field func- 
tional (9) up to the collapse point aN^^^ /auo = —0.61 (full 
line); the dashed line is the linear prediction (8). 

In order to take into account the effects of the inter- 
actions, we consider the mean field description based on 
the Landau- Vlasov equation 



^ + v^-^.^-i^— .^ = (11) 
dt dr dr dv 2 m dr d\ ' 



where /(r, v, t) is the distribution fimction, n = J fd^v is 
the atomic density. Eq.(ll) describes the dynamics of a 
normal weakly interacting gas in the coUisionless regime. 

An approximate solution of Eq.(ll), is obtained by 
making a scaling ansatz for the distribution function 



/(r,v,i) = /o(r(t),v(t)), 



(12) 



where /o is the equilibrium distribution, fi(t) = ri/hi 
and Vi{t) = biVi — biVi. Under the scaling assumption 
the velocity field u{r,t) = J vfd^v/n takes the simple 
form Ui — biVi/bi. This ansatz has been recently used by 
Guery-Odelin [11] to investigate the effect of the interac- 
tion on the collective oscillation of a classical gas in the 
coUisionless regime. 

The equations for the scaling parameters bi can be ob- 
tained by multiplying (11) by fj and Vi and integrating 
in phase space. Making use of the equilibrium properties 
of the distribution function, after some straigthforward 
algebra one finds 



bi + Jibi - + -xw- 



1 



0, (13) 
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where x is the ratio (4) evaluated at equiUbrium. The 
second term in (13) describes the restoring force of the 
external oscillator potential, the third one originates from 
the kinetic energy, while the last term, linear in x, ac- 
counts for the effects of mean field interaction. 

An immediate application of Eq.(13) concerns the 
study of the oscillations of the gas. By linearizing the 
equations around equilibrium (6j = 1) one finds, in the 
presence of isotropic harmonic trapping (wj, =Wz = wq), 
the result: 

WM = 2a;oVl + 3x/8, wg = 2u;o^l - 3x/4, (14) 

for the frequencies of the monopole and quadrupole os- 
cillations which coincide with the results already derived 
in [10] using a sum-rule approach. 
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FIG. 2. Aspect ratio as a function of time for the expan- 
sion of the normal (lower curves) and superfluid phase (upper 
curves) for A = 0.1 and two difi^erent values of the parameter 
X: X = (full liiie) and x = —0.4 (dashed line). 

The equations describing the expansion are obtained 
by suddenly removing the second term of Eq.(13), orig- 
inating from the trapping potential. In the study of the 
expansion we are interested in the case of anisotropic 
trapping. In particular we will consider the case of cigar 
shaped traps. For high deformations (A = w^/wj, ^ 1) 
Eq.(13) yields the asymptotic result 



6^-a;,^(l-^x)^^ 



showing that the aspect ratio 
i?±(t) ^ 1 
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approaches a value smaller than 1 if the interaction is 
attractive (x < 0). In Eq.(17), i?_L and Rz are the radii 
where the atomic density vanishes (Thomas- Fermi radii). 



The results of the numeric integration of the equations 
of motion (13) are reported in Fig. 2 and 3 as a function 
of time for the choices x = and x = —0.4. 

We address now the problem of the expansion of a su- 
perfluid Fermi gas. As already anticipated we will make 
use of the hydrodynamic equations of superfluids. Those 
equations have been already used to describe the collec- 
tive oscillations of a superluid trapped Fermi gas [12] in- 
cluding its rotational behaviour [13]. The hydrodynamic 
equations arc applicable if the healing length is much 
smaller than the size of the sample, which implies that 
the energy gap should be larger than the oscillator ener- 
gies , . This non trivial condition implies that the 
whole system behaves like superfluid. Furthermore the 
hydrodynamic equations are applicable up to excitation 
energies of the order of the energy gap. In the problem 
of the expansion it is crucial that the system remains 
superfluid in the first instants when the hydrodynamic 
forces provide the relevant acceleration to the expanding 
atoms. One expects that this condition be satisfied if the 
initial temperature is small enough. 




FIG. 3. Aspect ratio as a function of time for the expan- 
sion of the normal (lower curves) and superfluid phase (upper 
curves) for A = 0.3 and different values of the parameters X- 
X = (full line) and x = —0.4 (dashed line). 

The hydrodynamic description is based on the equa- 
tion of continuity 
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and on the Euler equation 



m^u + V (^^He{n) + Vfco(r) + ^mu^^ = (19) 

where niein) is the chemical potential of a uniform gas 
calculated at the density n and u is the velocity field. 

If the equation of state is a power law (/x^e oc rf) these 
equations admit the simple scaling solution 
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and the Thomas-Fermi radii evolve according to the law 

Ri{t) = Ri{0)bi{t). In this case, it is immediate to show 
that, during the expansion, the scaling parameters obey 
the coupled differential equations 
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for highly elongated configurations (A <^ 1). For 7 = 

1 (corresponding to a Bose-Einstcin condensed gas) the 
equation for the radial motion is integrable analytically 
and one finds the result b_i_{t) = (1 + ^1*2)^2 [5]. 

To describe the expansion of superfluid Fermi gas we 
will use the same equation of state (2) as for the normal 
phase. The case of a very dilute gas is also described by 
a power law with 7 = 2/3 (first term in (2)). For A = 0.1 
and 0.3 the solution is given by the full upper line in Fig. 2 
and 3 respectively which show that the deformation of 
the trap is inverted in time and the aspect ratio R^/R^ 
reaches asymptotically a value significantly larger than 
1 [14]. Superfluidity has hence the effect of distributing 
the release energy in a strongly asymmetric way along the 
axial and radial directions. It is worth noticing that the 
same scaling equations (22), with 7 = 2/3, are obtained 
for a classic gas in the collisional regime [15]. 

In the more general case (2), a useful approximation 
to the solution of the hydrodynamic equations, based on 
the scaling ansatz (20,21), is obtained by multiplying the 
Euler's equation (19) by rin{r) and integrating over the 
spatial coordinates. Using the equation of state (2), one 
finally obtains the following set of differential equations 
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with X defined by Eq.(4). Equations (24) differ from the 
anlalogous equations (13) holding in the normal phase. 
By linearizing Eqs.(24) around bi = 1 one gets, in the 
case of a spherical trap, the result uiq = V^loq for the 
quadrupole frequency [12], which, contrary to (14), is 
independent of the interaction term in 

The predictions of Eqs.(24) for the expansion of the 
gas are reported in Figs. 2, 3 and show that the inclusion 
of the interaction term in the equation of state affects the 
expansion of the superfiuid only in a minor way. 



In conclusion we have shown that the expansion of a 
superfiuid Fermi gas, being governed by the equations of 
hydrodynamics, differs in a crucial way from the one of 
a normal gas in the coUisionless regime. From a theoret- 
ical point of view several questions remain to be investi- 
gated: among them, the effect of large scattering lengths 
[16] on the equation of state and the role of collisions 
which, under certain conditions, might give rise to a hy- 
drodynamic regime, and hence to anisotropic expansion, 
also in the normal phase. Finally one should develop the 
formalism at finite temperature where both the normal 
and superfluid components are present. The resulting 
bimodal structure in the expanding cloud is expected to 
be affected by the transfer of atoms from the superfiuid 
to the normal component during the first stage of the 
expansion. 
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